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ABSTRACT 

In order to investigate the origin of the interstellar turbulence, detailed observations in the CO J = 
1-0 and 3-2 lines have been carried out in an interacting region of a molecular cloud with an H II region. 
As a result, several 1,000 to 10,000 AU scale cloudlets with small velocity dispersion are detected, 
whose systemic velocities have a relatively large scatter of a few km s _1 . It is suggested that the cloud 
is composed of small-scale dense and cold structures and their overlapping effect makes it appear to 
be a turbulent entity as a whole. This picture strongly supports the two-phase model of turbulent 
medium driven by thermal instability proposed previously. On the surface of the present cloud, the 
turbulence is likely to be driven by thermal instability following ionization shock compression and 
UV irradiation. Those small scale structures with line width of ~ 0.6 km s _1 have a relatively high 
CO line ratio of J =3-2 to 1-0, 1 < R 3 -2/i-o ^ 2. The large velocity gradient analysis implies that 
the 0.6 km s _1 width component cloudlets have an average density of 10 3-4 cm -3 , which is relatively 
high at cloud edges, but their masses are only < 0.05 M . 

Subject headings: Turbulence - ISM: structure - ISM: clouds - ISM: kinematics and dynamics - stars: 
formation 



1. INTRODUCTION 

Understanding the physical mechanism of contraction 
and condensation of molecular cloud cores toward star 
formation is a long standing issue. To date, two ma- 
jor paradigms have been proposed. First is that cores 
are in a magnetically subcritical state, and as the mag- 
netic field dissipates by ambipolar diffusion they co llapse 
quasi-statically (magnetic model; Shu et al. 1987). The 
other model is that the interstellar medium (ISM) is dom- 
inated by supersonic turbulence and the collision of the 
turbulent flow forms temporary density enhancements, 
some o f which undergo gravitati onal collapse (turbulent 
model; |Mac Low & Klessen|2004|). The magnetic model is 
favored by some theoretical studies because the magnetic 
field can be modeled to have arbitrary strength and con- 
figurations (aligned or disordered). It is, however, very 
difficult to be investigated observationally, being a 3- 
dimensional vector field. On the other hand, turbulence 
is ubiquitous in interstellar space. Almost all molecular 
line emission observed in molecular clouds have larger 
line widths than those expected from their thermal sound 
speed. A difficulty of the turbulent model is its unknown 
nature, namely the origin of the turbulence. The inter- 
stellar turbulence must have a much shorter dissipation 
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timescale than lifetime of dense cores because its super- 
sonic flow forms shocked regions and the energy is rapidly 
lost by radiation. The turbulent model requires a driving 
mechanism to keep itself over a long timescale. 

These two models are supported by some observational 
evidence. Low density diffuse clouds in general have fil- 
amentary structures, where aligned magnetic field lines 
(parallel or perpendicular to the long a xis) are often ob- 
served (e.g., Mye rs fc Goodman||199l| . These indicate 
that the ISM is frozen onto the magnetic field lines rather 
than bound by self- gravity. Molecular cloud cores are, on 
the other hand, reported to be dominated by turbulence, 
and tho se forming stars tend to have relatively small line 



widths ([Tachihara et al. 2QQ2[ ). This suggests that as 
interstellar turbulence dissipates, cores tend to get grav- 



inter stellar turbulence dissipates, cores tend to get grav- 
itationally bound, condensed, and form stars. It is also 
shown that the amount of turbulence is different among 
regions possibly due to difference in the environment; for 
example, cores in Taurus have much smaller line widths 
than those in Ophiuchus North and Lupus where star 
formation is not active and the times cale of core evo- 
lution is sup posed to be longer (e.g., |Tachihara et al. 



2000a 2002). Hence, turbulence is suggested to be one of 
the key elements controlling star formation activity. For 
these reasons, investigating the origin and driving mech- 
anism of turbulence is vitally important for determining 
the timescale and necessary conditions of star formation. 
A clue may be in the environment, namely, many OB 
stars are distributed in the Ophiuchus North and Lu- 
pus regions where cores have relatively large line width. 
A strong UV radiation field is known to affec t kinemat- 
ics and morp hology of molecular clouds (e.g., Lefioch & 
Lazareff|p94l ). 

A promising theoretical model for the origin of in- 
terstellar turbulen ce was proposed by Koyama & Inut- 
suka (2000, 2002) and following studies (Yamada et al. 
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20071 H noiie fc Inu tsuka 2008, 2009a; Hennebelle fc Au- 
dit||2007| |Hennebelle et al.||2QQ7| ) Their model (here- 
after the two-phase medium model) suggests that shock 
compressed layers of the interstellar diffuse gas undergo 
thermal instability and fragment into small scale (< 1000 
AU) structures composed of denser cold neutral medium 
(CNM). They are supposed to coalesce to form larger 
structures (~ 10000 AU). These cloudlets are embed- 
ded in warm neutral m edium (WNM), and have ra ndom 
motions (see Fig. 1 of |Koyama fc Inutsuka 2002| ). The 
velocity is supersonic with respect to the sound speed 
of the CNM (T ~ 50 K), bu t subsonic with respect to 
that of WNM (T ~ 8000 K) ( [Field et al.|[l969| |Koyama 
fc Inutsu ka 2000). Because of the two-phase medium, 
turbulence is sustained for relatively long timescales. 

In order to corroborate the idea of the two-phase 
medium model, we intend to detect small scale structures 
of the molecular cloud and investigate how interstellar 
turbulence is driven from the clouds' morphological and 
kinematical structures and their physical properties. We 
make detailed observations of CO J =1-0 and 3-2 at the 
surface of the molecular cloud LDN 204 facing an H II 
region where ionization shock compression is expected to 
be taking place. 

2. TARGET REGION LDN 204 

LDN 204 is a filamentary cloud complex facing the Sh 
2-27 H II region (r ~ 5 pc), that is excited by the near- 
est O star, C Oph (SpT = 09.5V, d = 140 pc). This 
cloud complex was entirely surveyed by the NANTEN 
telescope in 12 CO (HPBW = 2.7), and it was found 
that t he cloud is accelerated by illuminating UV radi- 
ation ([^chihar^F^ [2^ The 
cloud complex horbours embedded dense 3ores discov- 
ered in C 18 0, while star formation is n ot active and no 
associa ted young stars have been found (Tachihara et al. 
2000a). 



We carried out CO J =1-0 observations at 115 GHz 
with the NRO 45 m telescope (HPBW = 15", corre- 
sponding to 2000 AU). The interface region of the cloud 
surface of 11' x 22' area centered at (R.A., Dec.)j2ooo 
= (16:46:45.0, -12:21:30.0) is surveyed with the BEARS 
multi-beam receiver employing the On-The-Fly (OTF) 
mapping mode. We repeated the vertical and horizontal 
scanning many times and summed the data afterward 
in order to avoid scanning effects. We also employed 
scaling for correcting the beam efficiency by multiplying 
the beam-specific factors provided by the Nobeyama ob- 
servatory. The telescope pointing was checked every 2 
hours by obtaining 5-point observations toward nearby 
SiO maser sources. 

Throughout the observing runs, the wind speed has 
been below 15 m s _1 . About one fifth of the observations 
were taken under slightly windy conditions as wind speed 
was > 10 m s _1 . We thus divide the data into two to see 
how the spatial resolution gets worse under the windy 
conditions. No significant difference between windy and 
calm conditions can be seen, and eliminating the data 
with windy conditions worsens the S/N ratio of the final 
image. We therefore decided to use all the data. 

The data taken with the OTF mapping is reduced 
on the NRO data re duction software called "NOSTAR" 
( Sawada et al.|[2008 ). The baseline is subtracted by fit- 
ting to a linear function. The data is re-gridded into a 
7'/ 5 grid spacing by convolving with the first order Bessel 
function x the Gaussian function, resulting in the final 
data cube with an effective spatial resolution of 19'/ 5. For 
each spectrum of the 7'/ 5 grid cell, a final signal sensitiv- 
ity of T rms = 0.37 K with a velocity resolution of 0.1 km 



is achieved. 



3.2. ASTE observations 



The line width of C 18 J =1-0 observed by the NAN- 
TEN telescope in this region is typically as large as 
~ 0.74 km s _1 , significantly larger than that in Tau- 



rus (0.49 km s A m average; 
strong UV radiation is suggested to have compressed the 
molecular cloud leading it to form dense cores, while the 
streaming motions of the low density gas imply consid- 
era ble kin etic energy input to the ISM QTachihara et 
aL]|2000b|). Noteworthy is that ( Oph is a run-away 



Toward three selected regions where peculiar small 
scale structures are detected, we made follow-up CO 
J =3-2 observations with the Atacama Submillimeter 
Telescope Experiment (ASTE). It has a heterodyne re- 
ceiver working at 345 GHz with typical noise tempera- 
Onishi etTal Tl996b" The ture °^ ^sys ~ 400 K (double-sideband) at the zenith. It 



star traveling through the region ejected from the cen- 
ter of the Upper Sco subgroup of the Sco OB association. 
Therefore the entire cloud complex is under the influence 
of the UV radiation from £ Oph within a few x 10 5 years. 

Another interesting aspect is that the filamentary 
cloud is penetrated by magnetic fiel d lines perpendic- 
ular t o the long axis of the filament (McCutcheon et al. 
1986). The cloud complex is suggested to be controlled 
by magnetic fields, while no further investigation toward 
the dense region has been reported. 

Because of the conditions mentioned above, in partic- 
ular its proximity and the strong UV radiation, we chose 
the interface region of LDN 204 and Sh 2-27 as the best 
target for the present study. 

3. OBSERVATIONS 
3.1. ^5 m telescope observations 



is equipped with the 1024-channel digital backend that 
has a frequency resolution of 125 kHz, which corresponds 
to the velocity resolution of 0.11 km s _1 . Its 10 m dish 
enables us to have comparable spatial resolution to the 
45 m telescope with the HPBW of 22". The two tar- 
get regions are centered at (16:46:45.0, —12:22:50) and 
(16:46:36.5, -12:28:30) with sizes of 4' x 4' and 3' x 3', 
respectively, and are mapped with the OTF mode. The 
pointing error was corrected every 2 hours measuring the 
offset from a strong point-like source in CO. The calibra- 
tion for the beam efficiency is done by observing spectra 
of M17SW (R.A., Dec.)j 20 oo = (18:20:23.1, -16:11:43) 
assumin g that the peak r adiation temperature (T^) is 



69.6 K ([Wang et al.||1994| . The data are reduced with 
NOSTAR, similar to the 45 m data, resulting in the final 
data cube with the effective spatial resolution of 24" 2. 
The rms noise temperature (T rms ) of the spectra of each 
7" 5 grid cell are 0.13 K after calibration. 

4. RESULTS 

4.1. Spatial and velocity structures in J =1-0 

Fig. [I] is the peak Tmb map of the CO J =1-0 line 
in the cloud boundary region of LDN 204. The west 
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side of the cloud is facing the H II region Sh 2-27 and 
the cloud surface is illuminated by UV flux from ( Oph 
located in the direction of north west. The edge of the 
cloud does not appear to be a smooth surface, but has a 
complex morphology, for example, there is a prominent 
clumpy structure protruding to the west in the middle of 
the observed area. The velocity field is also remarkably 
complex as seen in the channel maps of Fig. [2j 

Particularly a few prominent features can be recog- 
nized such as clumpy, pillar-like, and arc-like features as 
denoted by A, B, and C in Fig. |2j respectively. The 
clumpy structure is especially noticeable in the peak 
Tmb map as slightly detached from the main cloud. 
The clump (A) has a radius of ~ 45", correspond- 
ing to r ^ 6300 AU. Using the empirical relation of 
X = N(R2)/W(CO ) = 1.56 x 10 20 cm" 2 (K km s" 1 )" 1 
( Hunter et al.|1997| ), the total mass of the clump is esti- 
mated to be 0.05 M®, where N(H.2) is the column density 
of molecular hydrogen and W(CO) is the integrated in- 
tensity of the CO J = 1-0 line. Similar morphological 
structures have been discovered at interaction regio ns of 
the cloud surfac es with H II regions like M16 (e.g., Hes- 
ter et al. 1996[ ). They are referred to as "Evaporated 
Gaseous Globules" (EGGs). The size of the clump is, 
however, much larger than those of EGGs that range 
from 150 AU to 1800 AU. The pillar-like feature (B) 
is visible in the channels around Vlsr ~ 5.4 km s _1 
and overlapping with other velocity components. It is 
roughly pointing to the direction of the UV radiation 
source of ( Oph. The arc-like feature (C) is protruded 
from the main cloud by ~ 7000 AU, and if it is assumed 
to be a bent cylinder, its radius is r ~ 3000 AU. 

4.2. Velocity dispersions of the small scale structures 

The CO J =1-0 line profiles at the peaks of the 3 
prominent features and a typical line profile in the mid- 
dle of the cloud are investigated. Fig(3] illustrates that 
the lines at the 3 positions are fitted to single or double 
gaussian profiles and the central velocities and velocity 
dispersions are estimated. For positions A and C, the 
line emissions are fitted to single gaussians whose Full 
Width Half Maximum (FWHM) line widths are 0.57 km 
s _1 and 0.51 km s _1 , respectively. At position B, the 
line has clear double components, which have FWHM of 
0.63 km s _1 and 0.59 km s _1 . The central velocities are 
clearly different ranging from 3.8 km s _1 to 5.3 km s _1 . 
On the other hand, in the middle of the cloud, the spec- 
tral line shapes are not well fitted to the gaussian profile, 
but more irregular. They seem to be saturated at typi- 
cal peak antenna temperature of ~ 6 K, and consist of 
multiple velocity components overlapping on the line of 
sight forming large velocity dispersion. The composite 
average spectrum over the entire observed area is illus- 
trated by the black line in Fig. 3] The gaussian fitting to 
it gives a much larger line width of ~ 1.5 km s _1 than 
those of the small-scale structures (hereafter 0.6 km s _1 
component). 

The observed line width is composed of thermal and 
non-thermal components expressed as AV^ S = AVj.^ + 
AVf| T , where the thermal component of the velocity 
dispersion at temperature T is expressed as AT^h = 
^8 In 2 • kT/m where k is the Boltzman's constant and 
m is the mass of the molecule. For the cases of T = 



20,60,150 K, AVth of the 12 CO molecule is calculated 
to be 0.2, 0.35, 0.55 km s _1 , respectively. Because we 
do not expect temperatures as high as 1000 K for molec- 
ular clouds, typical line widths of ~ 1.5 km s _1 have 
been interpreted as turbulent rather than thermal mo- 
tion. The velocity dispersion of ~ 0.6 km s _1 implies, on 
the contrary, that the cloud small structures are nearly 
thermalized with temperature of ~ 150 K, or they have 
relatively small non- thermal components. 

These results imply that the molecular cloud is not 
made of a uniform material, but composed of small-scale 
cloudlets whose internal velocity dispersion is as small as 
the thermal motion of CNM, while the relative motion 
of the individual cloudlets is much larger. Because many 
small-scale cloudlets are overlapping on the same line of 
sight in the main cloud, the spectral lines obtained in the 
middle of the cloud are widened and exhibit larger line 
width mimicking supersonic motion of the entire mate- 
rial in the cloud. The relative motion of the cloudlets is 
in fact in the order of 1-2 km s _1 , which is supersonic 
with respect to the temperature of CNM (T ~ 50 K), 
but subsonic with respect to WNM (T - 8000 K). This 
interpretation is c onsistent with the two-phas e medium 
turbulent model of jKoyama fc Inutsuka (2002). 



4.3. CO J =3-2 results 

The data of ASTE CO J =3-2 and 45 m J =1-0 
are compared in detail for the purpose of investigating 
the gas density and temperature of the small-scale cloud 
structures. In order to make a point-to-point compari- 
son, the 45 m data are spatially convolved with a gaus- 
sian beam of 22" HPBW so that the 2 datasets have the 
same spatial resolution. 

Fig. [4] illustrates the peak Tmb ratio of R 3 _ 2 /i-o for 
the 2 regions where (a) the clumpy and (b) the arc-like 
features are detected. In region a, the clumpy feature 
entirely has #3-2/1-0 > 1 an d at the peak position al, 
-R3-2/1-0 ~ 1-4. The west side of the clump facing the 
H II region has even higher ratio of i?3_2/i-o ~ 2 at posi- 
tion a2. The west side of the clump has a steep intensity 
gradient and is illuminated by strong UV radiation, so 
the gas temperature is expected to be higher than the 
eastern side. Both of the J =1-0 and 3-2 lines at posi- 
tions al and a2 have line widths of ~ 0.6 km s _1 , while 
in the periphery at a3 and a4, their line widths are sig- 
nificantly larger than ~ 1.5 km s _1 (Fig. pa). In region 
b, the arc-like feature shows a similar trend of i?3_2/i-o 
having higher values in the west side. At the peak po- 
sition bl, # 3 _2/i-o is close to unity. The J =3-2 lines' 
intensity and shapes resemble the J =1-0 lines in general 
as seen at positions bl, b2, and b3. 

5. DISCUSSION 

5.1. Physical properties of the small scale molecular 

cloud 

In the observed region, there are small scale structures 
detected near the main cloud boundary, i.e., the 0.6 km 
s _1 width component whose typical size is 6000 AU < 



/ < 12600 AU. According to the theory of Koyama 



Inutsuka ( 2002[ ), thermal instability of shock-compressed 
W1N1V1 (T^SDOO K) forms thermalized cloudlets of a few 
100 AU made of CNM (T - 50 K) that are embedded 
in the WNM. They have random motion with a velocity 
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Figure 1. Pseudo color image of peak Tmb map in CO J : 
regions where CO J =3-2 data are taken with ASTE. 

dispersion of a few km s _1 and coalesce with each other 
to form larger structures up to several xlOOO AU. Our 
results are consistent with this scenario from the mor- 
phological and kinematical point of view. In order to 
further corroborate the theory, estimation of gas density 
and temperature are attempted in the following. 

The 0.6 km s _1 width component is well traced by both 
of the CO J =1-0 and 3-2 lines. The gas density and 
temperature can be constrained from the intensities of 
multiple lines by the large velocity gradient (LVG) anal- 
ysis, which requires density, temperature, column den- 
sity, abundance, and velocity gradient. For the calcula- 
tion, we adopt flowing assumptions setting the density 
(n(H2)) and the kinetic temperature (T^ n ) as free pa- 
rameters. 1) H2 density (n(H2)), CO abundance ratio 
to H2 (Z), and column density (N (CO)) have the rela- 



1-0 observed with the NRO 45 m telescope. The two red squares denote the 



tion of N(CO) = Zn(H 2 )/ where I is the path length of 
the cloud along the line of sight. 2) The path lengths 
(I) are the same as the apparent diameters of the clouds 
(2 x r), i.e., spherical or axial symmetry. 3) The velocity 
gradients are the same as the CO J =1-0 line widths di- 
vided by the path lengths (AV/l). 4) The CO abundance 
is constant as Z ~ 3 x 10 -5 for the entire clouds. The 
last assumption is, however, unlikely because the present 
region is under strong UV radiation and a considerable 
amount of CO molecules is dissociated to be C I, par- 
ticularly at the cloud surface. We therefore tested the 
calculation altering the CO abundance with 3 different 
Z values of 5 x 10~ 5 , 1 x 10~ 5 , and 5 x lO" 6 . As a 
result, it is found that the response of R 3 _ 2 /i-o to the 
same n(H.2) and T^ n changes with Z by a factor of 3 or 
less, while the velocity integrated CO intensity, W(CO), 
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Figure 2. Channel map of the CO J = 1-0 emission for the same region as Fig. ^ Each map shows CO intensity integrated over 0.3 km 
s velocity width as shown at the top right of each map. Three prominent structures (the clump, pillar and arc) are designated by A, B, 
and C, respectively, in addition to a middle part of the cloud by D. The data used to create this figure are available in the online journal. 



changes up to an order of magnitude. 

Figs. [6] illustrate the results of the LVG calculations 
showing J R 3 _ 2 /i_o as a function of n(H.2) and T k i n by gray 
scale and contours. Fig. [6^i is for the case of AV = 0.6 
km s _1 and I = 12600 AU, corresponding to the clumpy 
structure, which demonstrates the line ratio is in the 
range of 1 < R 3 _ 2 /i-o < 2 for the entire region, while 
the integrated intensity W(CO) in the clumpy structure 
is between 2.5 and 4.0 K km s _1 for the J =1-0 line. 
The parameter set of n(H 2 ) and T^ n that meets these 
observed conditions does not exist in this figure with 
T kin < 300 K for the case of Z = 5 x 10~ 5 . For the 
lower abundance cases, these conditions are satisfied with 
the parameter in the shaded areas of Fig. [6^i. For the 



case of Z = 1 x 10 5 , the density and kinetic temper- 
atures are restricted to be 1000 < n(H.2) < 3000 cm -3 
and 100 < T kin < 300 K; and for Z = 5 x 10~ 6 , they 
are 2500 < n(H 2 ) < 6000 cm" 3 and 40 < T kin < 120 
K. These are warmer and denser compared to typical 
dark clouds detected in 12 CO. From the size and thus 
constrained density, the mass of the clump is given as 
0.008 < M < 0.05 M , which is below the brown-dwarf 
mass limit. Since these masses are well below the Jeans 
mass at their current densities, these cloudlets are not 
expected to collapse gravitationally. Their further evolu- 
tion may include mutual coalescence, condensation from 
surrounding medium, "evaporation" (from cold gas to 
warm gas), and photo-dissociation. 
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- average AV = 1 .48 km/s 
-AAV =0.57 km/s 
-B AV= 0.63, 0.59 km/s 
-C AV=0.51 km/s 




of the 0.6 km s 1 component is relatively high compared 
to typi cal dark clouds, but they are still much colder than 



Figure 3. Spectra of the CO J = 1-0 line at the 3 peak positions 
of the clumpy (A in blue), pillar-like (B in green), arc-like (C in 
magenta) structures, and the middle of the cloud (D in cyan) as 
shown in Fig. [2] The black curve illustrates the average CO J = 
1-0 spectrum of the entire observed region. The dotted curves 
are the gaussian profiles fitted to those spectra, whose line width 
(FWHM) are indicated in the top-left of the figure. 

Fig. [6)3 is for I = 6000 AU and AV = 0.6 km s" 1 as for 
the arc-like structure, illustrating the similar trend, but 
using the same n(H 2 ) and 7ki n , #3-2/1-0 is calculated 
to be higher and W(CO) is lower than Fig. [6^, yielding 
lower kinetic temperature for the same R 3 _ 2 /i-o value 
and require higher density for the same IV (CO). Note 
that in the area on the left side of the turn-over point 
of the i?3_2/i-o contours, both lines are optically thin. 
As seen in Fig. [4]o, the arc-like structure exhibits higher 
-R3-2/1-0 than unity for the entire structure, and the 
peak ratio is ~ 2 around the position b4. The inte- 
grated intensity of the J =1-0 line in this structure is 
in the range between 2.0 and 3.5 K km s -1 . Hence, as 
for Fig. [6^l, there is no region meeting the observed con- 
straints with Tkin < 300 K and Z — 5x 10 -5 , while we ob- 
tained a satisfying region for the lower abundance cases. 
For the case of Z = 1 x 10 -5 , the density and kinetic tem- 
peratures are constrained as 2000 < n(H 2 ) < 5000 cm -3 
and 50 < T kin < 150 K; and for the case of Z = 5 x 10~ 6 , 
4000 < n(H 2 ) < 11000 cm" 3 and 25 < T kin < 60 K. 
These imply that the arc-like structure is cooler and 
denser than the clumpy one. From the derived den- 
sity range, the mass of the arc-like structure is also con- 
strained to be very low-mass as 0.005 < M < 0.03 M . 

If we use the empirical relation of A/"(H 2 ) = X • W(CO) 
with an assumption of X = 1.56 x 10 20 cm -2 (K km 



s -1 ) -1 ( Hunter et al.||l997 ), the average density is thus 
derived only from W(CO) and I as ~ 3300 cm -3 and 
~ 6000 cm -3 for the clumpy and arc-like structures, re- 
spectively. Under the circumstance with strong UV ra- 
diation, smaller Z and thus larger X are expected than 
above, yielding larger density, and therefore consistent 
density estimations are plausible. Nevertheless, these es- 
timates are uncertain under such an extreme condition. 

Although the mean kinetic temperature of molecular 
gas along the line of sight towa rd C Oph was estimated 
to be 54 K from H 2 absorption ( Savage et al.|1977 Liszt 
et al.|[2009| ), small scale structures may have large tem- 
perature variations. The estimated kinetic temperature 



WNM. [Tachihara et al.| (|2000bJ) reported that the entire 
LDN 204 cloud complex is a 20-pc long and 2-pc wide fil- 
ament with an estimated mass of 1100 M , yielding the 
average density of ~ 250 cm -3 . The above calculations 
reveal the clumpy and arc-like structures to be 4 to 40 
times denser than the cloud average. 

5.2. Other observational evidence of small-scale cloud 
structures 

From observational points of view, some studied have 
been published atte mpting t o investigat e small-scale 
structures of ISM. |Sakamoto| ([2002b 



2) and |Sakamoto fc| 
11-scale structures at 



Sunada (2003) detected similar small- 
the edge oThigh-latitude clouds (MBM 54 and MBM 
55) and Heiles Cloud 2 in Taurus, respectively, although 
their observations were by one-dimensional strip scan. 
The typical one-dimensional size is ~ 10000 AU, compa- 
rable to our 0.6 km s _1 structures, while their FWHM 
line widths are ~ 1.8 and ~ 2.8 km s" 1 for MBM 54 and 
MBM 55, and c± 2 km s _1 for HCL2, respectively, signif- 
icantly larger than the present results. Without strong 
UV radiation, these regions are supposed to have smaller 
pressure and remarkably low density of n(H 2 ) < 100 
cm -3 is estimated, although no detailed estimation of the 
physical cloud parameter was performed. These clouds 
with different velocity widths may have been caused by 
d ifferent initial condition s . 

|Falgarone et al.| p009 ) investigated even smaller-scale 
structures using the IKAM-PdB interferometer in the Po- 
laris Flare, and detected elongated structures of ~ 300 
AU thickness. They tend to be seen in the CO J =1-0 
line wings with line width of 0.1-0.4 km s _1 . Because 
some of the elongated structures form parallel pairs with 
different velocity, it was suggested that the small struc- 
tures are thin layers of CO cloud with velocity shears 
driven by turbulent flow. 

At high galactic latitude, small-scale CO structures 
have been discovered by singl e- dish survey and by inter- 
ferometry. |Heithausen| ( |2002[ ) serendipitously discovered 
small-scale (< 1') structures with AV ~ 0.8 km s _1 with 
the IRAM 30m telescope. The clouds are resolved into 
smaller-scale (~ a few 10") structures with AV ~ 0.4 km 



s^ 1 b y the PdBI interferometric observation ( Heithausen 
2006). Assuming the distance of 100 pc, their sizes and 
densities are < a few xlOO AU and > 20,000 cm -3 , re- 



spectively. By a follow-up large-scale survey, Heithausen 
(2006) claimed that such small clouds are common in 



the ISM, but their total mass amounts to only < 1% of 
the whole gas mass. The size and AV of the small-scale 
structures ranging from ~ 100 AU to 10 4 AU includ- 
ing an ensemble cloud have a correlation of AV oc r 0,3 , 
where r is the effective radius of the cloud. Our sample 
of the 0.3 km s _1 and the 0.6 km s _1 width components 
r oughly fit to the corr elation. 

Ingalls et ah] ( [2007] ) observed MBM-12 with OVRO 
interferometer and detected 1-5 milli-parsec scale clouds 
in the CO line wings. They discussed that the velocity 
field of the cloud is made of an ensemble of tiny diluted 
turbulent cells. 

In the course of the natal cloud sur vey of the young 
stellar cluster TW Hya association, Tachihara et al. 
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Figure 4. Pseudo color image of the peak Tmb ratio of R 3 _ 2 /i-o around the clump (a), the arc (b), respectively. Overlaid are contours of 
peak T a J = 1-0 from 1.0 K (= 2.7a) with 0.4 K steps. The coordinates of the center positions are (16:46:45.0, -12:22:50) and (16:46:36.5, 
-12:28:30), respectively, and the sizes of the maps are A' x 4/ and 3 X x 3 X , respectively, as denoted by the red boxes in Fig.^ 
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Figure 5. Spectra of CO J = 1-0 (black) and 3-2 (red) at the positions marked in Fig. pi The CO J = 1-0 data are convolved with 
a gaussian beam of 22" HPBW so that the 2 data cubes have the same spatial resolution, j^he intensity and the spectral shapes of the 
emission lines of AV ~ 0.6 km s _1 are similar for both transitions except for al and a2. The data used to create Figures [l] and [5] are 
available in the online journal. 

As for the physical nature of the ubiquitous inter- 
stellar turbulence from diffuse to dense gas, the Kol- 
mogorov model has been supposed to be a plausible 
mechanism. According to the model, turbulence is driven 
by a larger-scale wave such as a supernova explosion, 
stellar wind, and molecular outflows, and cascades into 
smaller-scale structures. This cascade results in an en- 
ergy spectrum expressed by a power-law function as 
E(k) oc /c -5 / 3 where E is the energy of the turbulent 
flows and k is the wave number . Some observational re - 
sults support this model (e.g., Armstrong et al. 1995), 
however, the driving mechanism is known to have a dif- 
ficulty that the large-scale supersonic flows cause shock 
dissipation and most of the input energy is rapidly lost by 



(I2QQ9D detected small clouds by the NANTEN telescope. 
They have complex velocity structures and patchy mor- 
phology, barely resolved by the low spatial resolution ob- 
servation (2. / 7), but the line width is relatively narrow as 
~ 1 km s _1 . Some of them are, on the other hand, re- 
vealed to have small distance (< 100 pc), indicated by 
the optical Na absorption lines. This ensures that the 
total molecular mass derived from CO within the sur- 
veyed area of 1.9 deg 2 is < 3 M@. These results suggest 
that such small-scale cloud structures are common in the 
ISM. 



5.3. Nature of the small-scale structures and origin of 
interstellar turbulence 



8 



Tachihara et al. 
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Figure 6. Intensity ratio diagrams of CO R 3 _ 2 /i-o as a function of molecular hydrogen density (n(H2)) and kinetic temperature (T^[ n ) 

shown by the gray scale and contours assuming the CO abundance of 5 x 10 -5 . Contours are from 0.2 to 2.8 with steps of 0.2, and the 
bold contours indicate the line ratio of unity and 2. The red and green lines are contours of 0.2, 1.0 and 2.0 for the cases of lower CO 
abundance with Z = 1 X 10 -5 and 5 x 10 -6 , respectively. Used parameters are (a) I = 12600 AU and AV = 0.6 km s _1 and (b) I = 6000 
AU and AV = 0.6 km s _1 . The calculated integrated intensity of CO J =1-0 under the given condition are shown by dotted contours of 
2.5 and 4.0 K km s _1 for (a) and 2.0 and 3.5 K km s _1 for (b), respectively. The shaded areas indicate the parameter ranges meeting the 
observed values (see text). 



radiation. In order to solve this problem, many magneto- 
hydrodynamic (MHD) simulations have been attempted 
with a help of the magnetic Alfven wave. Nevertheless, 
it is reported that the MHD turbulence decays within 
a timescale of the order of the flow crossing time even 
with strong magnetic fields (|Mac Low et al. 1998 Stone 



|et al.|1998| ), because the cascade of the MHD turbulence 



happens within an eddy turnover time ( Cho & Lazarian 
2005]). 

According to the two-phase medium model, on the 
other hand, the molecular clouds in turbulent interstellar 
medium are built up from small scale structures formed 
by thermal instability and their coalescence. The ISM 
consist of two-phase gas, WNM and CNW, that can 
weaken the shock dissipation of the supers onic motion 



an d prol ong the decay time of turbulence (Inutsuka et 
al.||2005|). In the p Ophiuc hi cloud, CO and C I gas "Is 



reported to be well mixed ( Kamegai et al.||2003 ), which 



supports the idea of the two-phase medium. II a cloud is 
interacting with an H II region like the present case, the 
cloud surface is exposed to the ionization shock, that pro- 
duces turbulence via the thermal instability. The detec- 
tion of the small scale structures with cold temperatures 
and high densities is strong support for the two-phase 
medium model. 
In addition to t he results of density d i stribu ti on of the 



origin a l wor k by Ko yama fc Inutsuka| (2002), Yamada 
et al.| (|2007|) further investigated the two-phase model 



including the radiative transfer calculation to estimate 
the temperature distribution as well as calculation of the 
line ratios of CO and C II emissions. They suggested that 
the small-scale CNM are formed from WNM via thermal 
instability having a small dense structure with minimum 
temperature of ~ 50 K, and the estimated R 3 _ 2 /i-o is 



between 2 and 3, slightly higher than our result. Their 
condition is, however, set to be lower density (n(H 2 ) < 
3000 cm -3 ), and for the case of higher density and lower 

temperature, R 3 _ 2 /i-o is expected to b e smaller. 

No te that the numerical simulation s of |Koyama fc InutT 



suka| ( 2002 ) and Yamada et al. ( 2007[) begin wit h uniform 
WNlvr ^ompressed by passage of a shock wave. |Koyama 



Inutsuka 



(2002) showed that if the pre-shock gas is 
essentially one-phase diffuse WNM, the thermal insta- 
bility induced by shock compression results in forming 
two-phase turbulent medium. However, the present case 
of our observed molecular cloud (CNM) is supposed to be 
clumpy (multi-phase) medium that preexists even before 
the UV ionization takes place. This is because the UV 
source of £ Oph is a run-away star, it has been irradiat- 
ing the cloud only for a few xlO 5 yrs, wh ich is shorter 



th an the timescale of thermal instability (Tachihara et 
aLl|2000b ). In addition, the molecular CO gas seems to 



be well mi xed with atomic C I g as even in dense parts of 
the cloud (iKamegai et al.|[2003| , the ISM in this region 
is supposed to be clumpy two-phase medium even before 
the interaction wit h C Oph. O n the other ha nd, using 
MHD simulations, |Inoue et al] ( |2009b[ |2011| ) revealed, 
however, that passage of shock wave through the two- 
phase medium drives the ISM turbulence owing to vor- 
tex creation by shock-cloud interaction. This mechanism 
does not require thermal instability, and hence, the tur- 
bulence is driven more effectively with shorter timescale 
than the gas cooling time, and even shorter than the 
crossing time of ( Oph. It is also anticipated that un- 
der the strong UV radiation the ISM should have larger 
turbul ence due to larger input energy causing larger line 
width ( jTachihara et al.||2000a| |2002| [Gritschneder et al.| 
20091. 



Origin of Turbulence in Molecular Clouds 



9 



In order to further study the origin of interstellar tur- 
bulence, high resolution interferometric observations in 
multi-transition CO lines at cloud boundaries are re- 
quired. Also for the investigation of physical properties 
of WNM, extensive C I observation in the submm wave- 
length is highly desired. Observational results brought 
by ALMA combined with ACA would enormously im- 
prove our understanding in the near future. 

6. SUMMARY 

In order to investigate the origin of interstellar turbu- 
lence and verify the theoretically proposed thermal insta- 
bility model of two-phase medium, we have carried out 
CO J =1-0 and 3-2 observations toward a boundary be- 
tween LDN 204 dark cloud and an H II region Sh 2-27. 
The main conclusions are as follows. 



5. From the LVG analysis, it is found that the 0.6 
km s _1 width components have average density of 
~ 10 3-4 cm -3 , much higher than the entire cloud 
average density. Owing to their small sizes, masses 
of the 0.6 km s _1 width components are estimated 
to be 0.005 < M < 0.05 M , well in the domain of 
brown-dwarf mass. 

6. These results strongly support the theoretical ther- 
mal instabili ty model of two-phase medium first 
proposed by |Koyama fc Inutsuka| ( |2000[ |2002| ) . We 
conclude that the interstellar turbulence is likely to 
be driven by thermal instability and the small-scale 
structures are the building blocks of the turbulent 
molecular clouds. 



1. By detailed investigation of spatial and velocity dis- 
tribution of CO J =1-0, small-scale cloud struc- 
tures have been discovered. 

2. They have characteristic morphologies such as 
clumpy, arc-like, and pillar-like structures, which 
appear in some velocity channels with typical size 
scale of a few xlOOO AU to 12000 AU and the ve- 
locity dispersion of AV ~ 0.6 km s _1 . The sys- 
temic velocities of these 0.6 km s _1 width compo- 
nent clouds have scatter of a few km s _1 . 

3. The spectral line profile taken in the middle of the 
main cloud, and the composite line profile of the 
entire cloud resemble single component broad-line 
spectra with velocity dispersion several times larger 
than the typical sound speed. These facts suggest 
that the cloud is composed of tiny cloudlets with 
small internal velocity dispersion. Because of their 
relative motion and overlapping effect on the line 
of sight, molecular clouds are observed to have su- 
personic turbulence as a whole. 

4. Follow-up observations in CO J =3-2 for these 
small structures revealed that the intensity ratio 
of i? 3 _2/i_o is close to unity for the 0.6 km s _1 
width component, which shows minimum kinetic 
temperature of T^ n > 25 - 100 K, relatively higher 
than typical dark clouds. 
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APPENDIX 

INTERSTELLAR DUST COMPONENT OF THE CLOUD 

Besides the molecular and ionized gas components, dust is an important content of the ISM. The emission, absorption 
and scattered light by the interstellar dust can provide morphological, geometrical and temperature information of the 
cloud. 

An Ha image was obtained toward the entire cloud complex of LDN 204 with VYSOS-6 A, a 150 mm telescopes 
of the Universitatssternwarte Bochum in Chile. Ten images each with an integration time of 2 minutes were taken 
during March 24th and 25th 2012. The images were taken with an Ha filter whose central wavelength and the FWHM 
are 656 nm and 5 nm, respectively, and in dither-mode for bad-pixels removal. After dark and bias subtraction and 
flatfielding the images were combined using a min/ max- reject ion. The reduced image was then binned (8 x 8) and 
low-pass filtered. In the observed region, the extinction of Ha is not as prominent as in the main filament of LDN 204, 
but visible as a dark nebular near the edge of the spherical H II region Sh 2-27 (Fig. [7]). This ensures that the entire 



cloud is located in front of the Stromgren sphere and gas is accelerated toward us by the UV light (Tachihara et al. 



2000b). 



The global distribution of the CO emission is compared with the dust thermal emission in 100 fim by Tachihara 



et al. (2000b). They are in good agreement in general, but the 100 fim emission is enhanced at the cloud surface 



irradiated by the UV light. The dust temperature distribution is represented by the color of far-IR thermal emission. 
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Figure 7. Ha image (2? 25 x 2? 25) of the LDN 204 cloud around the target region taken by the VYSOS-6 A telescopes at the Univer- 
sitatssternwarte Bochum in Chile. The regions observed by the 45m telescope and ASTE are designated by green and magenta rectangles, 
respectively. 

Fig. [8] is the 3-color composite image of the entire region of LDN 204 produced with the Improved Reprocessing of the 
IRAS Survey (IRIS) dataQ The cloud surface is heated by UV radiation and the dust has relatively high temperature 
as demonstrated by the enhanced 60 /im emission over 100 /im. Th e local high temperatur e causes the enhancement 
of 100 /im emissivity at the cloud surface relative to the CO line (Tachihara et al.||2000b|). The dense parts of the 



cloud are, on the contrary, in low-temperature indicated by relatively strong emission m 100 /im. 

Recent preliminary release of the Wide-field Infrared Survey Explorer (WISE) dats[^] is compared with the CO 
distribution. Fig. [9] shows 3-color composite image in the region around LDN 204 using the band 4 (22 jam) in red, 
band 3 (12 /im) in green, and band 2 (4.5 /im) in blue. The 22 /im emission is supposed to come from small dust 
particles not in thermal equilibrium, while the 1 2 /im one is rather domin ated by the emission from poly cyclic aromatic 
hydrocarbons (PAH) excited by UV light (e.g., |Compiegne et al. 2011). These mid-IR emission in general delineate 
the cloud surface facing the H II region, while tne emission in 12 /im appears to be more spre ad than in 22 /im. On 
the other hand, the dense part in the cloud traced by C 18 J =1-0 ( |Tachihara et al.||2000a[ ) appears to be a dark 
lane in mid-IR, unlike at 100 /im, slightly offset to the far side from ( Oph. The dust in the cloud is supposed to have 
high column density and low temperature in the dark lane. The CO J =1-0 distribution of the present study traces 
the edge of the cloud where the column density is not so high but mid-IR is bright in emission. Besides the cloud's 
edge, diffuse mid-IR emission is extended outside the region detected in CO. It suggests that considerable amount of 
molecular gas is photo- evaporated by the UV radiation but the warm dust grains remain at the cloud surface. This 



The 



IRIS 



data 



obtained 



from 



http://irsa.ipac.caltech.edu/data/IRIS/ 
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Figure 8. IRIS 3-color composite image of the LDN 204 cloud, 100' x 100' around the target region (12 /im in blue, 60 /im in green, and 
100 iim in red). The regions observed by the 45m telescope and ASTE are designated by green and magenta rectangles, respectively. 

feature is, on the contrary, not clearly visible in the absorption feature of Ha, implying that the warm dust cloud has 
relatively low column density. 
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